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Abstract 

We calculate the average elasticity considering non-diffractive and single diffrac- 
tive interactions and perform an analysis of the cosmic-ray flux by means of an 
analytical solution for the nucleonic diffusion equation. We show that the diffrac- 
tive contribution is important for the adequate description of the nucleonic and 
hadronic fluxes in the atmosphere. 



It is well known that the evolution of the nucleonic cosmic ray component is controlled 
by two physical quantities related to high energy hadron interactions: the interaction 
mean-free-path Xf^^, which is inversely proportional to the inelastic proton-air cross 
section, <7j^~ air , and the average elasticity (x) p_air , the fraction of energy retained by the 
incident particle after a collision. It was shown in Refs. Jl], that, when one supposes 
the interaction mean-free-path and the mean elasticity as energy-dependent quantitites, 
the analytical solution for the nucleonic diffusion equation in the atmosphere is given by 

F N (E,t) = N E-^ +1) exp 

where t is the atmospheric depth and NqE'^ 1 ^ is the primary differential spectrum. 

In previous papers |TJ, we have discussed the importance of the energy dependence 
of the leading particle spectrum, which was done through the mean inelasticity. As a first 
approach, we have neglected the diffractive interactions in the proton-air collisions. By 
diffractive interactions we mean processes like 

a + b — > a + X 

where particle b is excited to a system X with the same quantum numbers and character- 
ized by an invariant mass M. This kind of process is called single diffraction (for details 
see, for instance, Ref. In general, diffractive interactions are neglected in cosmic-ray 
physics (mostly in analytical calculation of cascades), because their contributions are a 
priori considered to be very small. 

In this paper, we consider that the leading particle distribution has two contributions, 
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namely non- diffract ive (ND) and single-diffractive (SD), and we show that, although the 
latter is really small, its effect is relevant for an accurate description of cosmic-ray fluxes. 

In spite of its simplicity, Eq. (1) has two parameters, (x) p ~ a,n and Af n _air , which must 
contain all dynamic aspects of the hadron collisions occurring in the atmosphere. There- 
fore, in the following we shall describe how these aspects can be taken into account and, 
mainly, how to compute diffractive effects in the leading particle distribution. 

In Eq. (1), the interaction mean-free-path is given by 

K^(E) = 2 p i:'°l (g/cm 2 ), (2) 
< (mb) 

with the p-a.ii inelastic cross section calculated here by means of the Glauber model 

<- mi = J d 2 b {1 - exp [-o&T(b))} , (3) 
where b is the impact parameter and T(b) is the nuclear thickness 

r+oo 

T(b) = / p(b,z)dz (4) 



given in terms of the nuclear distribution p(b,z) (see Ref. 0). For cr^t, we use in Eq. 
(3) the best fit of UA4/2 Collaboration §. 

Now we need to establish (a;) p_air in order to calculate the nucleonic flux. We begin 
by writting the partial average elasticities, SD and ND, as 

/ \ JO dx I \ JO x dx /r\ 

{x)sd = ; (^)nd = , (5) 

0"SD 0"ND 
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which are correctly normalized to compound the average elasticity in pp collisions by 



VP , p P _ PP N 
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In the above expression, we are assuming that af£ = o"sd + a NE>5 where a(£ is given 
by the Landshoff parametrization, of£ = 56s~ ' 56 + 18.16s 0,08 [[/J. 

In the whole calculation, whose results we shall show farther on, we have used two 
models from which we borrowed the ND and SD distributions. For the non-diffractive 
elasticity ((^)nd), we use the Interacting Gluon Model (IGM), revised by Duraes et al. 
with the purpose of including semi-hard interactions, responsible for mini-jets events. In 
Ref . H , we have shown that the correspondent non-diffractive leading particle distribution 
produces reasonable results for pp total cross section and for p-air inelastic cross section. 

In order to include the diffractive contribution, we use the leading particle distribution 
of the Covolan-Montanha model [0] , which reads 



do-sv f dVffl, ^ (?) 



dx J dtdM 2 
where, for pp interactions, the invariant cross section is given by 



rfVsp_ = (S(3 p G p (t)) 2 (M2) 
dtdM 2 16tt IPpV ; ' 



(8) 



°jp p {M 2 ) = 3/U(^(M 2 )>, (9) 
with (3 P = 2.502 GeV~\ £ = 0.0764 GeV, G p (t) is the electric form factor of the proton, 
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ccjp(0 = 1-08 + 0.25£, where t is the squared four-momentum transfer, and (r 2 (M 2 )) = 
12.75 + 0.84 In M 2 . We remind the reader that, in diffractive processes, we have x = 
1 - M 2 /s. 

As we intend to calculate hadronic interactions in the atmosphere, (x) pp must be 



corrected to include air effect. This is done by the procedure given in Ref. [[10 1 



fin 



( X )P-™ = J2 Pn({x) P Ti (10) 



n=l 



where 



p. = J ^ (6) di: 



p— air 



and 



P n (b) = -i K t T(6)]"exp [-o&T(&)] . (12) 



Here P n is the probability of n-fold collisions of the primary nucleon inside the nucleus 
and n max , the maximum number of collisions, is roughly given by 2.3A 1 / 3 ||10|| . 

The diffractive contribution included in Eq. (6) comes from inelastic interactions 
among the incident hadron and the nucleons inside the struck nucleus. In addition to 
this contribution, it is necessary to consider the diffractive dissociation of the nucleus as 
a whole. For this reason, we add to < x >^~ air (once again, in a weighted way) a second 
component coming from nuclear diffractive processes, so that our final expression is 



p— air / p— air p— air\ 

<*> P - alr = ^Ir-^lD^ + [(Jm p'J SD } (XY-^. (13) 



In order to calculate (x)g D dir and ctsd^' we use an extension of the Covolan-Montanha 
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model to nuclear diffractive interactions which is made by a radial scaling f9|, i.e. multi- 
plying Eq. (|) by 



(14) 



(7-2)2 

where (r^(A))^ = 1.096 —0.41 corresponds to the atomic radius as a function 



of the atomic mass A, and (r^)? = 0. 197^12.75 + 0.84 Ins is the hadronic radius of the 
proton, both given in fermis. 

As experimental data for the nucleonic flux cover a large range of energy, 1 < E (GeV) < 
10 3 , we need to stablish a cutoff for our elasticity at laboratory energy E = 50 GeV; for 
lower energies, the (x) pp is kept constant and the (2) so is turned off. This is necessary 
basically because of two reasons: 1) for energies below 50 GeV, the elasticity given by 
the IGM drops steeply to very low values; 2) at these low energies the Covolan-Montanha 
model, which is based on the Triple- Pomeron formalism, is out of its validity range and 
does not hold. In order to avoid a discontinuity of the calculated flux due to this abrupt 
change in the elasticity, we use its value at 50 GeV as a constraint and calculate the cor- 
respondent (x) pp . Then, this value ((x) pp = 0.6125) is kept constant for lower energies, 
only being corrected to include air effect. 

In the Fig. 1, the calculated nucleonic flux is compared with experimental data mea- 
sured at sea level [|ll|, |12| , using the Ryan primary spectrum |13[], N = 2(cm 2 .s.sr.GeV)~ 1 



and 7 = 1.5. With only ND contribution, our solution underestimates the nucleonic flux, 
as shown by the dashed line in Fig. 1. The final result, including both contributions 
(ND and SD) and free of parameters, is shown by the solid line. One can see a significant 
improvement of our theoretical calculation, when diffractive effects are taken into account. 
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Recently, measurements of hadronic fluxes at sea level, in the energy range 1 < 
E (GeV) < 10 5 , were performed by the hadronic calorimeter of the KASCADE exper- 
iment |14]]. In order to compare the hadronic flux calculated by means of the formalism 



developed in this paper with those data, we borrow the KASCADE parametrization for 
the pion to nucleon ratio 



R = - +?r = 0.04 + 0.271og(£/GeV), (15) 
n + p 



and we correct with this factor our nucleonic flux to obtain the hadronic flux 

F H (E,t) = (1 + 2R)F N (E,t). (16) 

In the Fig. 2, the analytical solution of the hadronic flux is compared with KAS- 
CADE's experimental data, showing a good agreement between them. 

In conclusion, the diffractive contribution to the mean elasticity is taken into account 
by two mechanisms. The first one comes from the inelastic interactions among the in- 
cident hadron and the nucleons of the target, while the second considers the diffractive 
dissociation of the nucleus as a whole. The net effect of these two corrections is to increase 
the value of the mean elasticity. This enables us to perform a satisfactory description of 
the nucleonic and hadronic fluxes at sea level, without further parameters. 

We would like to thank the Brazilian governmental agencies CAPES, CNPq and 
FAPESP for financial support. 
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Figure Captions 



FIGURE 1: Nucleonic flux at sea level. Dashed line: only ND contribution. Solid line: 
ND and SD contributions. Open circles, Ref. [11]. Full circles, Ref. [12]. 

FIGURE 2: Hadronic flux at sea level. The solid line is our calculated hadronic flux. 
Stars, Ref. [14]. Open diamonds, Ref. [15]. 
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